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Abstract 

The electromagnetic productions of J/ip associated with light hadrons(LH) and leptonic 
pairs(^ + /i~, t + t~) at B factories are studied. We find that the direct electromagnetic produc- 
tion cross sections of J/ip(ip(2S)) associated with light hadrons is about 0.10(0.04) pb. The di- 
rect production cross sections of J /ip{ip(2S)) associated with is about 0.056(0.020) pb. If 
we include the contributions from ip(2S) decay, we can get the prompt cross section a[e + e~ — > 
J/ip + n+fJ,'] = (0.068 ± 0.002) pb, about (16 ± 5)% of the Belle data a[e + e~ -> J/i/j + X non _ c5 ] = 
(0.43 db 0.09 ± 0.09) pb, meanwhile the e + e~ —> J/tp + r + r _ process only contributes 2%. The 
prompt cross section a[e + e~ — )• J/ip + Light Hadrons] = (0.121 ± 0.005) pb is about (28 ± 8)% of 
the Belle data. 

PACS numbers: 13.66.Bc, 12.38.Bx, 14.40.Pq 
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I. INTRODUCTION 



Heavy quarkonium physics is an important ground to test quantum chromodynamics 
(QCD) both perturbatively and non-perturbatively. Non-relativistic quantum chromody- 
namics (NRQCD) factorization approach l| has achieved a series of successes in describing 
heavy quarkonium production and annihilation decay. However, there are still some pre- 
dictions which are less satisfactory. For more details, see a concise review Q. Among the 
problematic comparisons with experiment, the large discrepancy between theoretical predic- 
tions and experimental data on the charmonium production in e + e~ annihilation has drawn 
much attention recently. 

In recent years, the B factories has provided systematic measurements on charmonium 
production. Some results are puzzling, because of the large gap between the measurements 
and the theoretical predictions. For example, the cross section cr(e + + e~ —> J/if> + r) c ), mea- 



sured by Belle Collaboration [3] and Babar Collaboration Q], is almost one order-of-magnitude 
larger than the leading-order(LO) predictions {5- 7]. By introducing the QCD perturbative 



correction [8|, |9(, and in combination with relativistic correction lOi, lllfl. this discrepancy was 
largely resolved. Besides the challenges in the exclusive process, the large ratio of J /if) pro- 
duction associated with charmed hadrons is also confusing, which was measured by Belle j^] 

o[e + e~ — >• J /if) + cc 



a e T e 



J /if) + X] 



0.59ig;g ± 0.12. 



(1) 



In contrast to the leading-order NRQCD predictions, this ratio is only about 0. 1 [12Ml4| . The 
next-leading-order (NLO) QCD corrections were also introduced in J/if> inclusive production 
to resolve the discrepancy between experimental measurements and LO calculations. The 
NLO QCD corrections to e + e~ — > J/ibcc process enhance the cross section with a K factor 

_ nn 

16] . and only about 20 percent for e + e — > J /if>gg process [17J, [18j. So the 



of about 1.8 
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discrepancy is greatly alleviated. To check what role does the color-octet proce ss p lay, the 
NLO QCD corrections to color-octet J /if) inclusive production was calculated [19[. Com- 
bining the relativistic corrections to e + e~ — > J/ipgg[20\, it may imply that the values of 
color-octet matrix elements are much smaller than the expected ones which are estimated 
by using the naive velocity scaling rules. 

Another interesting topic is the e + e _ — > J /if> + X non _ c5 production. Most recently, the 
prompt J /if) production in association with charmed and non-charmed final particles was 



measured 21] 



a(e + e- J/ip + X)=(1.17 ± 0.02 ± 0.07)p6, 
a(e + e- J/ip + cc>=(0.74 ± 0.08^;^)p6 
a(e + e- J /if) + X non _ c5 )=(0A3 ± 0.09 ± 0.09)p6. (2) 



These processes are investigated in Ref. 20|, |22|, the results show that including both the 
0(a s ) radiative correction and the 0(v 2 ) relativistic correction, the color-singlet contribution 
to e + e~ — > J/ipgg has saturated the latest observed cross section e + e~ —> J/tj) + X non _ cZ 
measured by Belle. 



Aside from the above QCD process, the pure QED process should also be considered [23]. 
Especially, in our paper, we calculate virtual-photon- associated production a(e + e" — > 
J I4>1*) x -^(t* ~^ ll(or Light Hadrons)). The rest of the paper is organized as follows. In 
Section II, we will give the formulations of e + e~ — > J/ip + fi + fi~. In Section III , the QED 
production of e + e~ — > J/ij) + LH is discussed. In section IV, we will give the numerical 
results and discussion. Finally we summarize our results in section V. 



II. THE FORMULATIONS OF e+e J/tp + (j,+fj,- 

In NRQCD factorization scheme, the cross section of e + e~ 4 J/^ + /U + pT can be de- 
scribed as follows 

^l(e + (A; 1 )e-(A; 2 ) -> J/if>{2p) + n + (Pi) + ^M) 

= y^Yl EE^ Si; ^ 2 1 SS z)( LL ^ SS * I JJz)(^k | 1) x 

L Z S Z siS2 jk 

A(e+e- cf{p) + c s k 2 (p) + ^( Pl ) + iT{ P2 )) (3) 

where (3j;3k | 1) = 1/a/3, (si; s 2 | SS Z ), (LL Z ;SS Z \ JJ Z ) are respectively the color-SU(3), 
spin-SU(2), and angular momentum Clebsch-Gordan coefficients for cc pairs projecting out 
appropriate bound states. A(e + e~ —> Cj{p) + Ck{p) + H + (pi) + H'ipz)) is the scattering 
amplitude for cc production. The coefficient C$ can be related to the radial wave function 
of the bound state and reads 

C s = ±- | R S (0) | 2 • (4) 
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FIG. 1. The Feynman diagrams of e + e — > J/ipfi + [i 



We introduce the spin projection operators PssAPiQ) as 

Pss*(p,q) = ^2(si, s 2 \SS z )v(p - q; sjuip + q; s 2 ). (5) 



Expanding Pss z {p, q) in terms of the relative momentum q, we get the projection operators 
at leading term of q, which will be used in our calculation, as follows 

P 1Sz ( p ,0) = -Lf( Sz )(tf + M/2). (6) 

where M is the mass of the charmonium. It is two times of charm quark mass m in the 
non-relativistic approximation. The polarized cross section can be calculated by defining 
the longitudinal polarization vector as follows 

where 4p 2 = M 2 and n M = (1, —pj \ p |). 

The Feynman diagrams of e + e~ — > J/tpfi + fi~ are shown in Fig{TJ The process of final 
states with plus charge parity is depicted in diagram Fig{T](a, b) and denoted as C = +. 
The process of final states with minus charge parity is depicted in diagram FigJIJc, d) and 
denoted as C = — . The C = + process is dominant, and the C = — process is suppressed 
by a factor of 



here the logarithm term come from quasi-collinear divergence with m 2 ^ s. The s/M 2 term 
come the photon propagator. One can get / ~ 10~ 2 for e + e~ — > J l^p&yT at y/s = 10.6 GeV 
for B factories. 
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III. THE QED PRODUCTION OF e + e~ -»■ J/ip + LH 



Similar to the e + e~ — > J/ipfi + fi~ process, the fragment process is also dominant in 
e + e~ J/ip + Light Hadrons process. The fragment process can be considered as e + e _ — > 
7*7*, then the virtual photons fragment into J/ip and light hadrons respectively. By using 
the approach of the calculation of the hadronic part of the muon g — 2 241 ]. this process can 
be described as 

da QED [e + e~ -» J/ip + LH] da[e + e~ J/ip + y+if 



dm 2 LH 



dm 2 + _ 



x R 



had / 2 

rn l i+ l jr 



,(9) 



=m LH 



where 



R had^ 



a[e + e — > hadrons] 



(10) 



a[e+e ] m2 =A2 

because of the contribution of the C = — process is negligible, after subtracting the effect 
of 7* -)■ cc from R had (A 2 ) by a naive factor 4/36 (A - M cg _ T/ie ), we get 
da QED [e + e" ->• J/ip + LH] 



dm 2 LH 



da[e + e — > J /ip + 
dm 2 , _ 

MM 



X 



R had (m 2 



mm 



-0(m 



2 

M + M~ 



M. 



cc-The) 



where 9 is step function, M c5 _ The should be correspond to the cc threshold of Mj/^p, 
etc. The uncertainties should be discussed in the next section. 



(11) 
2M D , 



IV. NUMERICAL RESULT 



In numerical calculations, the parameters are selected as 25] : 

M M = 0.1057GeV, Mj^ = 3.0969GW, = W.QGeV, 

M T = 1.7768GeV, M fm = 3.686GeV, a = 1/132.33 (12) 



The table of R had (X 2 ) have been given in Ref . [24j| . We construct an interpolation of 
R had (A 2 ) corresponding to the table with first-order interpolation and setting R had (A 2 ) = 
when A is less than \ min in the table. 

The wave function at the origin can be extracted from the leptonic width T(V — > l + l~) 

2 

Tn 

|i?s(0)|2 = 4^ r[V ^ e+e1 - (13) 



The leptonic width of charmonium decays into e + e has been given in Ref. 25] 



r[J/V> -> e + e~] = 5.55 ± O.UkeV, 
T[ip(2S) -> e + e~] = 2.38 ± 0.04A;e\/. (14) 

When we calculate the prompt production cross sections of J/ip, we take into account 



the feeddown contribution from ip{2S) by B[tp{2S) J/ip + X] = (57.4 ± 0.9)%[25[ and 
ignore the contribution of the other charmonium. Then we can get the direct production 
cross section of J /ip(ip(2S)) associated with r + r~ and at B factories as 

a direct [e + e~ J/^j + ^ >~] = 56 ± 2 /6 

cr direrf [e + e- -> J/V> + t + t-] = 6.4 ± 0.2 fb (15) 

and 

^ectjg+g- ^S) + fi + ^} = 20 ± 1 fb 

a dtrect [e + e- -)• ^(2S) + r + T~] = 1.8 ± 0.1 /6 (16) 

Most of the uncertainties come from the uncertainty of leptonic width in Eq.( fl4|) . The 
others come from the effect of fine structure constant a and higher order QED corrections 
and so on. The QCD corrections have been taken into account in the leptonic width of 
J/V(V(2S)). 

The cross sections for C = — process is only 1.6%(1.0%) of that for C = + process 
in J/ip(ip(2S)) production associated with n + fi~. And the ratio is about 6.0%(3.9%) in 
J/ij}(ij}(2S)) production associated with t + t~. These results are in agreement with the 
estimation in Eq.flH]). So the contribution of C = — process can be ignored in the calculation 
of electromagnetic J/i/j(ijj(2S)) production associated with light hadrons. Finally, we get the 
direct production cross section of J/i/j(i/j(2S)) associated with light hadrons at B factories 
as 

°qed t e+e ~ -> J /*P + LH \ = 100 ± 5 fb 
^QED[e + e- if>(2S) + LH] = 36 ± 1 fb (17) 

here we choose M c5 _The = 2Mn in Eq.fflll. If we choose M c5 -The = Mj/^, there is a difference 
of —1/6. So the uncertainties from M c£ _ The can be ignored. Most of the uncertainties come 
from R had and the leptonic decay width. 




The energy distributions of direct J/ip(ip(2S)) production from the processes e + e _ — >■ 
J/ip(ip (2S)ll and e + e~ — > J /ip{%l){2S))LH are shown in FigfS]and FigJ31 respectively. Unfor- 
tunately, the endpoint peak was not measured in Ref . . The polarization of J/tp(ip(2S)) 
direct production for e + e~ — > V + fi + fi~ (V = J/ip, ip(2S)) process as a function of the 
energy of J/ip is shown in FigJU The polarization of J/ip production associated with light 
hadrons is similar to the results of e + e~ — > V + /U + /i~ process. The angular distributions of 
direct J/i[)(i/}(2S)) production from the processes e + e~ — > V + fi + fi~(V = J/ip,ip(2S)) and 
e + e~ y 27* — > V + LH {V = J/ip, ip(2S)) are shown in Figj5]and FigJHJ respectively. 

Now we give the prompt production cross sections of electromagnetic J/tp production 
associated with leptonic pairs and light hadrons as 



o%%d i e+e ~ ->• + + X] = 68 ± 2 /6 

[e + e- ->• J/tl) + t + t- + X} = 7.4 ± 0.1 fb 



prompt 
QED 



a%Z Pt l e+e ~ -> + LH) = 121 ± 5 fb. (18) 
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FIG. 3. The J ftp and i/j(2S) energy spectra of direct production processes e + e — > 27 — > V + 
LH (V = J/i/),ift(2S)). 
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FIG. 4. The helicities of J ftp and ip(2S) direct production processes e + e — > V + (V 
J/ip,ip(2S)) as a function of the energy of J/ijj and ip(2S). 



V. SUMMARY 

In summary, the electromagnetic productions of J/ip associated with light hadrons(LH) 
and leptonic pairs(/z + /z~, r + r~) at B factories are studied. We find that the direct elec- 
tromagnetic production cross sections of J /if) (if) (2 S)) associated with light hadrons is about 
0.10(0.04) pb. The direct production cross section of J /ip(ip(2S)) associated with /z 4 ' \i~ 
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FIG. 5. The angular distributions of direct production processes e + e — >• V + (V 
J/ifj,if>(2S)). Here is the angle between J/if)(ip(2S)) momentum and beam. 
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FIG. 6. The angular distributions of direct production processes e + e — > 27 — > V + LH (V 
J/ip,ijj(2S)). Here 6 is the angle between J/tp(ip(2S)) and beam. 



is about 0.056(0.020) pb. If we include the contribution from tp(2S) decay, we can get the 
prompt cross section a[e + e~ — > J / 'if) + /i" 1 " pT + X] = (68 ±2) pb, about (16±5)% of the Belle 
data a[e + e- -)■ J/4> + X non _ cZ ] = (0.43±0.09±0.09) pb, meanwhile the e+e" -> J/ip + T + T~ 
process only contributes 2%. The prompt cross section <j[e + e~ —tJ/if) + Light Hadrons] = 
(0.121 ± 5) fb is about (28 ±8)% of the Belle data. Unfortunately, the endpoint peak of 
energy distribution for J /if) electromagnetic production associated with leptonic pairs and 
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light hadrons was not measured in Ref. 21J. The polarization of J/ip electromagnetic pro- 
duction associated with light hadrons is transversal, while the polarization of Jfap inclusive 
production associated with light hadrons from QCD process is longitudinal 22j . We also 
notify that the charge parity of final states is plus for the QED process calculated in our 
paper. And it is minus for color singlet process e + e~ — > J/ip + gg and color octet process 
e + e~ —> J/ij) + g. 
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